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We investigate the transport through a nanoscale device consisting of a degenerate double-orbital 
Anderson dot coupled to two uncorrelated leads. We determine the thermoelectric transport prop- 
erties close to the one-electron regime and compare them to a corresponding single-orbital dot. The 
linear and nonlinear regimes are addressed, the latter via a non-equilibrium generalization of the 
non-crossing approximation based on the Keldysh formalism. Power output and efficiency in the 
Kondo regime are shown to be strongly enhanced through the presence of a second orbital. We 
predict an experimentally relevant optimal operating point which benefits from the concomitant 
increase of the Kondo temperature in the two-orbital setup. An approximation based on the trans- 
port coefficients and fulfilling the thermodynamic balance is proven to remain appropriate even far 
beyond the expected range of validity of such approaches. Finally, the double-orbital Kondo regime 
reveals itself as a promising candidate to avoid, at least partially, the generic dilemma between 
optimal thermoelectric efficiency on one hand, and fair power output on the other. 
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I. INTRODUCTION. 

Over the last decade, the thermoelectric effects in 
nanoscale devices have become a field of intense re- 
search^, driven by the desire to reconvert waste heat back 
into usable electric energy. In addition to the obvious 
benefits for applications in microelectronics, where in- 
creasing clock rates and circuit integration make waste 
heat evacuation and power consumption serious issues, 
the subject is also interesting from a purely fundamental 
point of view: due to the high, essentially quantum-dot- 
like confinement in the conduction channel, electronic in- 
teractions and many-body correlations are expected to 
govern the transport properties of nanoscale devices. The 
most prominent effect is the so-called Coulomb block- 
ade, which tends to deteriorate the electric conductance 
through the conduction channel. Conversely, the Kondo 
effect, arising from resonant spin- flip scattering between 
the conduction electrons in the leads and electrons con- 
fined in the dot-like channel, allows to overcome the 
Coulomb blockade at low temperatures, thereby restoring 
the device's electric conductance^^. 

The impact of Kondo physics on the thermo- 
electric pro pert ies has become a topic under active 
experimental^ and theoretical^^ investigation. We 
would like to highlight the beautiful experiments on 
single-walled carbon nanotubes in the presence of tran- 
sition metal impurities^ where a giant thermopower has 
been revealed; and on lateral quantum dots where Kondo 
correlations in combination with electron-hole asymme- 
try have been found to strongly enhance the thermopower 
signaP. 

Recent theoretical works^M] have been dedicated to 
thermopower generation by nanoscale devices, focusing 
on efficiency and power output outside the Kondo regime. 
Mahan and SofxP3] predicted earlier that large values of 



the thermoelectric figure of merit ZT require a transmis- 
sion function of vanishing width across the device which, 
in turn, implies a very weak coupling between localized 
channel state(s) and leads, and ultimately results in low 
conductance and power output. A similar trade-off was 
uncovered by Nakpathomkun et ali^ in low dimensional 
systems. 

In this paper, we advocate that the Kondo effect may 
alleviate this dilemma for a quantum-dot like conduction 
channel: its Abrikosov-Suhl-Kondo (ASK) resonance can 
be narrowed at will by rising the on-dot Coulomb in- 
teraction U, no matter how large the value of the lead- 
dot hybridization T is. This route seems viable provided 
that temperature and voltage-bias effects do not erode 
the underlying Kondo physics. Formerly considered as 
an almost exclusive low-temperature phenomenon, the 
last decade's experimental progress boosted the Kondo 
effect to temperature ranges as high as several hundreds 
of KelvirPS and produced the first functional nanodevices 
with Kondo temperatures around 30 K F 6 I 27 1 Improve- 
ments can also be expected from the degenerate two- 
orbital Kondo effect, where entangled orbital and spin 
degrees of freedom have been shown, experimentalbJ^H 
and theoretically^, to raise the Kondo temperature sig- 
nificantly with respect to the standard spin-| Kondo ef- 
fect. 

In the following, we calculate thermal and elec- 
tric currents through a doubly-degenerate quantum 
dot. The calculation is performed in a fully nonlinear 
manner within the framework of the Keldysh Green's 
function formalism. Specifically, we use an out-of- 
equilibrium ge neral ization of the non-crossing approxi- 
mation (NCApMil which enables the computation of 
thermoelectric efficiency and power output, in presence 
of finite voltage and temperature biases applied across 
the quantum dot. The benefits for the thermoelectric 
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transport properties in the two-orbital Kondo regime 
become obvious in comparison to the standard single- 
orbital Kondo effect. Finally, we suggest an experimen- 
tally relevant optimal operating point. 



II. MODEL AND METHOD. 

The device considered in the following consists of a cen- 
tral dot coupled to two uncorrelated leads. At simplest 
level, retaining only one or two relevant orbitals in the 
dot, the situation is readily described by an Anderson 
model, with a Hamiltonian given in standard notation 
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The first line concerns the single (to = 1) or doubly de- 
generate (to = 1,2) orbital with on-site Coulomb repul- 
sion U, the second line describes the left (L) and right 
(R) leads, and the last line accounts for the hybridization 
tunneling between dot and leads which conserves orbital 
and spin quantum numbers. In the double-orbital case 
the above Hamiltonian is 5i7(4)-degenerate; nevertheless 
our results do not hinge on this symmetry. 

The results are found to depend little on the specific 
shape and width of the lead densities of states (DOS), 
provided that they are smooth enough and that their 
widths are by far the largest energy scale in the problem. 
For the solution of the NCA equations, we may there- 
fore use a Gaussian spin-summed DOS, N a (e), charac- 
terized by the same half- width D for both leads. The 
fc-independence of the tunneling amplitudes yields a hy- 
bridization between dot and leads that is directly pro- 
portional to the respective lead DOS, T a (e) = TTt^N a (e), 
(where henceforth a symmetric coupling, Fr(e) = Ti,(e), 
will be assumed). The total hybridization strength at 
the mean Fermi energy, JL — \{jll + Mr): will serve as an 
energy unit: T = r L (/x) + T R (jL). 

As appropriate for the power generating regime, the 
temperature bias AT = Tl — Tr > across the device 
is applied in the opposite direction of the voltage bias 
— eVb = Ml — MR < 0. The nonlinear electrical and hot 
(left) lead thermal currents are given, respectively, by 

1 =\\ [fL(e)-Me)}r(e)de and 

Ql = Ml) lh(e) - Me)} r(e) de , (2) 

where h denotes Planck's constant and — e the elec- 
tronic charge. fa(e) = f( e ~ Ma) are the lead 



Fermi functions, and r(e) = jA(e)T(e), with A(e) — 
— 4lm[^ m a G m<y (e + iS)], the total dot spectral den- 
sity. The retarded Green's function G ma (e + iS) = 
((cma', c mo-)) ^ s obtained fro m a generalized Keldysh- 
based out-of-equilibrium NCA'^SIMJ. The NCA is known 
to reliably describe the low-temperature Kondo scale Tk, 
and to give accurate results for the dot spectral den- 
sit y dow n to temperatures of the order of a fraction of 
7 r | 33 | 34 |. The power output is given by the sum of the heat 
currents, P = Ql + Qr, which together with Eqs. ^ 
yields P = IVh- The thermoelectric converter efficiency 
rj = P/Ql is conveniently normalized in terms of the 
Carnot efficiency, t]q = 1 — Tr/Tl. In the following, we 
focus on power generation, i.e. positive /. 

The transport in the linear regime is obtained by a 
first order series expansion of the Fermi functions in 
Eqs. ([2]) characterized by the following transport coef- 
ficients: electrical conductance G(T), Seebeck coefficient 
S(T), and electronic thermal conductance K e (T), 



G(T) 
S(T) 

K C (T) 



1 h(T)-vI (T) 



eT 



1 

T 



h{T) 



h{T) 
I 2 i(T) 
Io(T) 



(3a) 
(3b) 

(3c) 



given in terms of the transport integrals /„ (defined for 
n = 0,1,2)2] 
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where r cq (e) = jA eq (e)T(e) has the role of an equilib- 
rium transmission function. Finally, retaining only the 
electronic contributions to the thermal conductance, the 
figure of merit is defined by ZT = S 2 GT/K C . 



III. RESULTS. 

Taking JL as the energy origin, the chosen parameters 
£0 = — 3.2r and U = 16T, are close to the one-electron 
Kondo regime, implying approximately half filling for a 
single-orbital (SO), and quarter filling for the degenerate 
double-orbital (DO) case. More precisely, in the explored 
temperature range, the average number of electrons on 
the dot increases with T, and lies in the interval [0.86 — 
0.9] for SO, and [0.76 - 0.9] for DO. 



A. Linear transport properties. 

First we discuss the linear thermoelectric transport 
properties, some of which confirm earlier results by 
Sakano et a/P*l in Fig. [l] the temperature dependences 
of G, S and ZT are plotted for the SO and DO casesP. 
In the low-temperature regime, the Kondo effect allows 
to overcome the Coulomb blockade due to the presence 
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FIG. 1: (a) Electric conductance G(T), (b) Seebeck coefficient 
S(T), and (c) figure of merit ZT vs. T/F, for e = -3.2r and 
U = 16r. Dashed lines refer to the SO and plain red lines 
to the DO cases. Vertical dotted lines indicate the respective 
Kondo temperatures. 



of a robust ASK resonance in the immediate vicinity 
of the Fermi level which provides the necessary spec- 
tral density and thus guarantees conductance. Defin- 
ing the Kondo temperature as the value where G(T) 
reaches half its maximum, we find w 0.023r for the 
SO case (with = 1 henceforth), and a much larger 

(2) 

T£> « 0.5r for the DO case. Such an enhancement is 
in line with e.g. Hewson's textbook!^ and should be rel- 
evant for potential applications. As obvious from the 
dashed curve in Fig. [lji, the low-temperature conduc- 
tance for the SO case saturates close to the spin-summed 
quantum of conductance 2e 2 /h, implying that the con- 
ductance contribution of each spin species is almost max- 
imal, as expected for a device close to half filling. This 
contrasts with the corresponding low-temperature con- 
ductance for the DO case (full red line): for T <C T^ , we 
observe a large plateau where the conductance is almost 
temperature- independent. The height of the plateau does 
not quite reach 2e 2 /h expected at quarter filling, but sat- 
urates at G ~ 1.2e 2 /h, as predicted for an occupancy 
(n-ma) ~ 0.19, by the Landauer formula in combination 
with the Friedel sum rule: G = 2e 2 /hJ2 m sin 2 (n(n mcr )) 
(see e.g. RefJ^j). The DO conductance plateau ends to its 



right with a soft bump whose origin follows directly from 
the behavior of the spectral density: the ASK resonance 
is found to lie slightly beyond the Fermi energy, such that 
thermal broadening effects first enhance the spectral den- 
sity at the Fermi level, while only larger temperature rises 
yield to the inevitable erosion of conductance. 

Fig. [TJd displays the thermoelectric power as a function 
of T, whose sign indicates the particle- or hole-like nature 
of the transport. Again, the Kondo effect is respon- 
sible for the decisive feature, i.e. the low-temperature 
minimum of the Seebeck coefficient: due to its narrow- 
ness and location right above the Fermi level, the ASK 
resonance is very efficient in capturing a fraction of the 
hottest electrons of the left electrode, thereby giving rise 
to the thermoelectric effect. Approaching Tk the Kondo 
resonance vanishes and the effect disappears. Upon a fur- 
ther temperature rise, contributions from occupied states 
dominate which results in a sign change in S. For even 
higher temperatures, the unoccupied states with ener- 
gies around Eq + U reverse the sign change (off range). 
As can be seen from Fig. [T}d, this negative extremum in 
the Seebeck coefficient is located at roughly a quarter of 
the Kondo temperature. With respect to the SO case, 
the DO minimum is not only shifted to higher temper- 
atures, but also benefits from a Kondo resonance with 
significantly enhanced spectral weight, yielding a mini- 
mum which is almost twice as deep as the SO minimum. 
Although the present interpretation is strictly speaking 
only appropriate for the linear regime, it is confirmed in 
its main lines by the fully nonlinear calculation presented 
later in this section. 

As a direct result of the aforementioned behavior of the 
Seebeck coefficient, the figure of merit, shown in Fig. [lj;, 
benefits from this situation: its maximum is strongly en- 
hanced in the DO case, with the maximum ZT sw 0.4 
occurring for T$. « 0.25T K 2) w 0.12 T, i.e well above 

the SO maximum ZT w 0.1 at T^\ w 0.006 T. Parts of 
this enhancement are also due to the dip occurring in the 
ratio K e /(GT) for temperatures close to Tk, indicating 
a violation of the Wiedemann-Franz law, as noticed ear- 
lier by Krawiec et alP®. Although the temperature for 
which the figure of merit is maximal yields roughly the 
optimal efficiency of thermoelectric devices, we want to 
stress that this criterion has to be verified beyond the 
linear framework. 



B. Nonlinear transport. 

The results of such a fully nonlinear calculation for 
the DO setup are illustrated in Fig. [2j displaying the 
transport properties as functions of Vb and AT, which 
is applied such that the cold (right) lead is at constant 
temperature, Tr T^ pt . As obvious from the figure, effi- 
ciency and power output both exhibit a maximum, albeit 
for different regimes: the maximal efficiency, closely re- 
lated to the ASK resonance, is reached for AT ~ Tr, 
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FIG. 2: Nonlinear transport for the DO case in the power 
generating regime for eo = —3.2 T, U = 16 T, and fixed 
Tr = 0.12 F as functions of voltage and temperature biases 
across the device, (a) Normalized efficiency ri/r/c- (b) Power 
output. The various lines correspond to approximations based 
on transport coefficients: the dot-dashed and dashed lines 
represent, respectively, the lines of zero power output and 
the ridge lines for efficiency and power output for the TCA; 
the straight dotted line corresponds to zero power output of 
prosaic linear response theory (see main text). 



while the power output is maximal for (probably ex- 
perimentally unachievable) large AT ~ 5Tr where the 
ASK resonance is strongly suppressed. Our calculations 
predict an optimal efficiency of 10% of the Carnot effi- 
ciency for Vb fa 0.05r/e and AT ps 0.1T. Although the 
power output is not yet maximal at optimal efficiency, it 
reaches appreciable values of the order of 2.13 10~ 3 T 2 /h 
which could be interesting for potential applications. For 
an order-of-magnitude estimate, we assume a hybridiza- 
tion r = 25 meV, yielding a power output of 52 pW at 
maximum efficiency which is reached for Vb = 1.2 mV, 
T R = 35 K and T L = 64 K. 

The latter results should be compared to a Coulomb- 
blockade regime of similar efficiency. Postulating the ex- 
istence of an adjustable Lorentzian peak in the spectral 
function and choosing its width and position such that 
the efficiency maximum matches our calculations - i.e. 
10% of the Carnot efficiency - we find that the power 
output in the DO Kondo regime surpasses that of the 
Coulomb blockade regime by about 30%. This increase 
constitutes a significant step towards the resolution of the 
dilemma between optimal efficiency and power. For anal- 
ogous reasons to those mentioned in the linear regime, the 
physical origin of the power increase in the Kondo regime 



is due to the narrowness and position of the ASK reso- 
nance. We stress that this feature is not eroded by the 
present moderate voltage biases^. This clearly contrasts 
with the Coulomb blockade regime, where the width of 
the resonance is directly proportional to the hybridiza- 
tion, resulting in the aforementioned dilemma between 
large T, necessary for large power output, and small T, 
required for efficiency. From the experimental point of 
view, another advantage is that the ASK resonance is au- 
tomatically anchored in the vicinity of the lead Fermi lev- 
els, and hence does not require any fine-tuning of the level 
position £q. As a result, the thermogenerator remains op- 
erational for £o varying from — 4.4r to 0.7r, i.e. a range 
three times wider than in the pure Coulomb blockade 
regime. Furthermore, in line with our expectations, the 
DO setup clearly outperforms the standard SO one. In 
the latter case, the maximum efficiency is at only 3% of 
the Carnot efficiency, occurring at Vb ~ 2.38 10 T/e 
and AT ps 9.6 10~ 3 T, and the SO power output does 
not exceed 6.62 10~ 6 T 2 /h. 

Finally, we want to address the question whether sim- 
ilar results could have been obtained within the frame- 
work of an approximation which solely uses the trans- 
port coefficients (TCA). These coefficients were calcu- 
lated in the previous section and depend only on the 
device's equilibrium properties. To this encP^, we cal- 
culate the associated power output p TCA = Vf,/ TCA = 
—GV 2 — SGVbAT, and heat currents Q T L GA = SGT L V b + 
K e {l + ZT L )AT~\[I TGA ] 2 /G and Q£ CA = ^ S GT R V b - 
K e (l+ZT R )AT-^[I TCA ] 2 /G. The transport coefficients 
in the latter expressions, S, G, K e and Z, are evaluated 23 
at the device operating temperature T = §(Tl + 2r). 
Note that each reservoir temperature is used for the heat 
flux leaving the corresponding electrode^. The above ap- 
proximation also accounts for the Joule heating of the de- 
vice (which is small in our case), and presents the virtue 
of fulfilling the thermodynamic balance P — Ql + Qr- 
A detailed comparison reveals that this approximation 
reproduces the results of the fully nonlinear approach, 
qualitatively on the whole plot range of Fig. |U while 
small quantitative discrepancies appear for AT > 4T R . 
This remarkable agreement is illustrated in the figure by 
overlaying, on the fully nonlinear results, the TCA lines 
of zero power output (dot-dashed) and the TCA ridge 
lines (dashed). The efficiency ridge line is giverP^ by 

Vl^/Vc = Wl + ZT - l)/(v/l + ZT + T R /T L ). At 
this point, it is worth mentioning that linear response 
theory implemented around one of the reservoir temper- 
atures - in this case Tr - yields-^ the straight dotted 
line of zero power output in Fig. [2] which rapidly departs 
from the fully nonlinear results. This situation lead us to 
the conclusion that the discrepancy between such a pro- 
saic linear response theory and the fully nonlinear results 
does not necessarily require a genuine out-of-cquilibrium 
approach - even if this is sometimes claimecP-^21 - an d, 
conversely, that carefully implemented approximations 
relying on transport coefficients should be able to repro- 
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duce the nonlinear results rather accurately for similar 
devices. 

IV. CONCLUSION. 

In this paper, we have studied the thermoelectric trans- 
port through a double-orbital quantum dot in the Kondo 
regime. We showed that transport properties and power 
generation are strongly enhanced by the interplay of or- 
bital degeneracy and Kondo physics. This setup outper- 
forms not only the usual single-orbital quantum dot, but 



also devices in the Coulomb blockade regime. An ex- 
perimentally relevant optimal operating point has been 
determined. It owes its stability to the Kondo effect and 
might even allow for a complete removal of the gate elec- 
trode. A comparison of the transport properties obtained 
from carefully implemented approximations relying solely 
on the transport coefficients with the fully nonlinear so- 
lution proves the robustness of the former up to large 
applied temperature bias for the systems under study, 
thus adding credibility to many works on nanoscale de- 
vices conducted within this framework. 
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